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Abstract 

We investigate the effect of macromolecular crowding on protein folding, using purely repulsive 
crowding particles and a self-organizing polymer model of protein folding. We find that the ther- 
modynamics of folding for typical a-, (3- and a//3-proteins are well described by an adaptation 
of the scaled particle theory (SPT). In this approach, the native state, transition state, and the 
unfolded protein are treated as effective hard spheres with radii approximately independent of the 
size and concentration of the crowders. The same model predicts the effect of crowding on the 
folding barrier and therefore refolding rates with no adjustable parameters. A simple extension 
of the SPT model, assuming additivity, can also describe the behavior of mixtures of crowding 
particles. 



In contrast with conventional laboratory experiments conducted under dilute conditions, 
protein folding in a cell occurs in a dense environment consisting of various other macro- 
molecules, commonly referred to as "crowders" [1] . Although the detailed interactions of the 
protein with the crowders may be quite complex, the primary physical effect of macromolec- 
ular crowding on the protein folding reaction is to reduce the volume available to the protein 



by that occupied by the crowders 



3|, \% |5|, [y, |7| . A complete theoretical understanding of 



the excluded volume effects will greatly enhance our ability to interpret experiments, as well 
as all-atom simulations, and to develop coarse-grained models [8|, of concentrated protein 
solutions. 

Several theories have been put forward to predict the effect of crowders, modeled as 
impenetrable repulsive particles, on the folding free energy. However, these theories provide 
strongly contrasting predictions for the quantitative and sometimes even qualitative effects 
of crowding. For example, by treating the folded and unfolded proteins as effective hard 
spheres, Minton utilized scaled particle theory (SPT) to estimate the change in folding free 
energy as the difference between the insertion free energy for the folded and the unfolded 
states. The SPT free energy of inserting a hard sphere of radius R in a hard-sphere fluid 
of particle radius R c is |g|, f5F = (3y + 3y 2 + y 3 )p + {9y 2 /2 + 3y 3 )p 2 + 3y 3 p 3 - ln(l - C ), 
where (3 = 1/ksT, T is the temperature, k-g is the Boltzmann's constant, y = R/R c , 
P — 0c/(l — 0c)j and C is the fluid volume fraction. This theory predicts a rather strong 
effect of macromolecular crowding on the folding free energy 10] , with a monotonic increase 
in stability with increasing crowder packing fraction <j) c . 

An alternative theory proposed by Zhou also uses SPT for the the effect of crowders on 



the folded protein, but the free energy of the unfolded protein is calculated using an e 



egant 



model of the unfolded chain as a random walk in the presence of a spherical trap [11] . For 
a Gaussian chain of radius of gyration R g , with z = R g / R c , the change in free energy is 
(3F = 3(fi c z 2 +6ir~ 1 ^ 2 (j) c z — ln(l— C ). This model predicts a much weaker effect of crowding on 
the unfolded state free energy because the unfolded polypeptide can access the void space 
between the crowders. In addition, a maximum in stability is predicted as the crowder 
packing fraction is increased. The extension of this theory to binary mixtures of different 
size crowders leads to an intriguing conclusion that there will be an optimum mixing ratio 



of the two components to achieve maximum protein stability 



12]. One may expect that 



different crowding theories will work better under certain conditions of crowder packing 



2 



TABLE I: Characteristics of the proteins considered. 



Protein a N [A] otj [A] a t [A] [A] 

prb 9.5 12.75 11.4 15.9 

proteinG 10.6 16.0 - 20.0 

TNfn3 13.0 21.5 - 28.6 



fraction and size. However, the boundaries for the validity of these theories in parameter 
space are relatively unknown. 

In this letter, we investigate the appropriate theoretical description of purely repulsive 
crowders, using molecular simulations of a coarse grained folding model in a bath of crow- 
ders. We consider two-state proteins from the three main protein structural classes: all-a 
(prb [13|), all-/? (TNfn3 [14J) and a/ (3 (protein G |15| k which are described using a self- 
organising polymer model (or "Go-like" model) Il7fl . In our model, each amino acid 



residue is represented by a single particle and a standard procedure is used to build the 
potential from the experimental native-state structure If]]. Interactions between the con- 
tacts present in the native state are treated as attractive and all others as repulsive, an 
approximation motivated by the funneled nature of the folding free energy landscape 19]. 
Previous studies have shown that this type of simplified model can indeed capture relevant 



features of protein folding |20j, such as mechanism [21] and kinetics 22, 23]. The repulsive 
interactions between a crowder and the protein or other crowders are given by the pair 
potential V(r) = e[a rc f/(r — a + cr re f)] 12 , where r is the distance between particle centers, 
e = 1 kcal/mol sets the energy scale, a is the hard core overlap distance, and <T re f = 6 A is a 
reference diameter: for a = a re f, V(r) reduces to a more familiar form. We define a between 
the pair (i, j) as, a = Ri + Rj, where Ri, Rj are the radii of either crowder particles R c , or of 
the various protein residues R p . We use Langevin dynamics simulations with a time step of 
10 fs and a friction coefficient of 0.2 ps _1 , using the BBK integrator [24j in the CHARMM 
simulation package [25j. Cubic periodic boundary conditions with a primary cell size of 100 
A were employed. To speed up equilibration at a given temperature, we use replica exchange 
moves every 30 ps between 12 replicas which are each biased using an umbrella potential of 
the form Vi(Q) = 0.5n(Q — Qi) 2 , where Q is the fraction of native contacts, < Qi < 1 is 
the target Q value for replica i and k = 300 kcal/mol is the force constant. The required 
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FIG. 1: Effect of crowding on folding free energy surface of protein G. The potential of mean force 
along the reaction coordinate Q, the fraction of native contacts, is shown for bulk conditions and 
with different crowder packing fractions <f> c for crowder size R c = 8 Aat 320 K. The curves have 
been shifted by an arbitrary constant to match the free energy of the folded state. Inset: Heat 
capacity C V (T). 



thermodynamic information is extracted from simulations with different umbrella potentials 



26. 



27]. 



and temperatures using the weighted histogram method (WHAM) 

For the model proteins considered here, Q is nearly an optimal coordinate for identifying 
transition states and also capturing the dynamics of protein folding with a diffusive Marko- 



vian model 



23 



28l . |29j. Here, we project folding in our systems onto Q, and use it to 



identify transition states as well as unfolded and folded free energy minima. To estimate 
the crowding- induced changes on the folding free energy surface, we construct the potential 
of mean force (PMF) as a function of Q, denned as (3F(Q) = - \n[P{Q)/AQ], where P(Q) 
is the equilibrium probability of observing configurations between Q and Q + AQ. Figured] 
shows representative free energy profiles in bulk (0 C = 0) and under crowded conditions 
(0 C > 0, R c = 8) for protein G. The destabilization of the unfolded state due to the presence 
of crowders is clearly visible from an upward shift in the curves near the unfolded basin 
(Q ~ 0.1). For high crowding packing fractions (0 C > 0.10), there is also a slight destabi- 
lization of the transition state (Q ~ 0.5) with respect to the folded basin (Q « 0.9). Fig. [1] 
(inset) shows the heat capacity C v curves for the same state points, as obtained from the 
WHAM analysis. The maximum in C v is defined as the folding temperature Tf. The shift 
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FIG. 2: Shift in protein stability with crowder size and packing fraction. Top panel: pair distri- 
bution functions g(r) between crowders (C) and the native (N) and unfolded (U) protein (C-C [U] 
is the C-C g(r) in the presence of U), for prb with R c = 8 A, <f> c = 0.05. Arrows indicate corre- 
sponding sums of hard sphere radii. Lower panels: change in folding free energy with crowding: 
upper: protein G, middle: prb, lower: TNfn3. Symbols: simulations (errors smaller than symbol 
size) for R c = 20 A (circle), 16 A (square), 12 A (triangle up), and 8 A (triangle right); dashed 
lines: predictions from Zhou theory [ll|; dot-dash lines: predictions from Minton theory with 
a\j = ^/5/3-Rg ; solid line: result of using a smaller hard sphere for the unfolded state, motivated 
by the form of the U-C g(r). 

in Tf toward higher temperatures with increasing <p c demonstrates the stabilizing effect of 
the crowders on the folded protein relative to unfolded. 

Figure [2] shows the effect of crowding on the free energy of folding, AAF = AFu-n(0c = 
0) - AFu_ N (0 c ), at T = 320 K for prb and protein G and T = 300 K for TNfn3, where 
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AFu_ N = -fc B Tln(/ Qt e-P p WdQ/ Jq^-^UQ), is the difference between the unfolded 
F\j and native free energy, Q$ is the location of the transition state along Q. We find 
a monotonic increase in stability with the packing fraction C for all proteins and crowder 
sizes. For a given <p c , the smaller crowders are more strongly stabilizing because there are 
effectively fewer voids of the size of the protein. The prediction of Zhou's theory using the 
folded protein radius and unfolded protein radius of gyration computed from the 
simulations (Table [I]) gives excellent agreement with the simulation data in the limit <p c —> 
(Figure HJ. For smaller crowders and higher packing fractions, however, the theory predicts 
too small a stabilization, most likely due to the neglect of excluded volume in the unfolded 
chain. 

On the other hand, the Minton theory predicts a stabilization much greater than found 
in the simulations over the full range of C . This theory treats the unfolded protein (with 
radius of gyration iff) as an "equivalent hard sphere" of uniform density, i.e. with radius 



y 5/3.Rg . However, as indicated by the pair distribution function g(r) between the unfolded 
protein and the crowders (Figure [2]), the unfolded protein is quite soft, so that a smaller 
choice of hard sphere radius may be more appropriate. Remarkably, we are able to fit the 
data over the full range of C and crowder size (which varies by a factor of 2) to the SPT 
model by using the radius of the unfolded state a\j as a single adjustable parameter. This 
radius is indicated on the plot of g(r) in Figure [2j At high C the smallest crowders are 
better able to penetrate the protein and the fit could be marginally improved by using a 
slightly smaller ajj. This treatment of the unfolded protein as a hard sphere is motivated by 
the highly successful Asakura-Oosawa theory [30j of polymer-colloid mixtures so that soft 
sphere interactions can be used between the polymer and crowder 3l| which in turn can 
be mapped onto a hard-sphere system. In our case, we find that the effective hard-sphere 
radius a\j of the unfolded chain is ~ 80% of R^. 

We have also studied mixed macromolecular crowding with binary A:B and ternary A:B:C 
mixtures for the a-helical protein prb. For A:B mixtures, keeping the packing fraction of 
the type "A" crowders (radius 8 A) fixed at <f>^ = 0.05, we varied the fraction of type 
"B" crowders (radius 12 A). The simulation results for the change in folding free energy 
for different mixtures are shown in Figure [31 To test whether the effect of mixed crowd- 
ing is simply additive and can be easily estimated from the pure crowder simulations, 
we also calculate the change in folding free energy from the following additivity ansatz, 
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FIG. 3: Scaled particle theory prediction for binary mixtures of 8 A and 12 A crowders. With 
the packing fraction of particle type A fixed to 4>^ = 0.05 and the fraction of particle type B, 
4>c = 0.05 - 0.20, the change in folding free energy of prb from simulations AAF s ; m is shown as a 
function of change in folding free energy AAF^d from the additivity ansatz (circle) as discussed in 
the text. The inset shows simulation data (circle) along with SPT model predictions (line) using 
the additive model. 

AAF add (0j, <pl, ...ific) — J^iXiAAFiiJ^nPl), where index i runs over iV different types of 
crowding particles, Xi = (j) l J J2j 4^ i s the fraction of crowder type i in the mixture, and (f) l c is 
the volume fraction of i. The predictions of this additive model AAF a dd for A:B mixtures 
are in extremely good agreement with the simulation results AA_F s ; m as shown in Fig. [31 
Moreover, for an A:B:C mixture of crowder radii 8, 12, 16 A with volume fraction 0* = 0.05 
of each component, the agreement between the simulation result (0.55 kcal/mol) and the 
additive model (0.54 kcal/mol) is remarkably good. Therefore, it should be possible to es- 
timate the effect of mixed crowding due to several types of repulsive crowder by utilizing 
single crowder results. Applying additivity to our SPT model predictions for single crowders 
also provides reasonable estimates for AAF (Figure [3] inset). Our results and the additive 
model for mixed crowding do not predict an optimal mixing ratio as expected from a pre- 
vious theory [ijj]. Indeed, for a mixture of crowders at given C , the greatest stabilization 
will occur when all the crowders are of the more stabilizing (smaller) type. 

From our explicit dynamical model of folding we are able to estimate folding rates di- 
rectly from mean first passage time calculations. Starting from at least 400 different initial 
coordinates drawn from an equilibrium unfolded ensemble (Q = 0.2) at a given packing 
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FIG. 4: Influence of crowding on folding kinetics. The effect of crowding on the folding barrier 
height of prb is estimated using k-gT ln(/cj? ulk / fcf) for crowders of radius 8 A (circle), 12 A (square) 
and 16 A (t riangle up); lines are the SPT predictions. Inset: Folding rates with binary crowder 
mixtures (A: 8 A radius, B: 12 A radius, <fif = 0.05). Results from direct simulations (circle) are 
compared with the predictions based on simulations with a single crowder size (square) and from 
SPT (line). 

fraction, we calculate the average time Tf = 1/kf taken to reach the folded state (Q > 0.9). 
We estimate the change in barrier height (AF|_u — F$ — F\j; F$ = F(Q$)) with crowding as 
ksT ln(fc^ ulk / fcf) , which is justified if the position of the folding barrier and the diffusion co- 
efficient along the reaction coordinate are unchanged (Figure H]). The folding barrier height 
decreases monotonically with increasing <p c : indeed, it is possible to predict the change in 
rates using SPT with no further adjustable parameters: we calculate the transition state 
radius oj directly from the Q umbrella simulations (Q = 0.55), and estimate the change in 
barrier height from SPT, using a$ in place of on- The effect of mixed macromolecular crowd- 
ing on kinetics can also be obtained from single crowder simulations by assuming additivity 
(Figure HI inset). 

In summary, we find that scaled particle theory provides an accurate description of the 
effect of macromolecular crowding on both folding stability and rates. The effect of purely 
repulsive crowders can be well-approximated over a wide range of crowder sizes and pack- 
ing fractions by treating the unfolded state as a hard sphere of fixed radius ajj. In all the 
cases considered, we find a\j fa 0.8Rg. Our results have several important consequences. 
For crowders of a single size, folding rate and stability will increase with increasing packing 
fraction monotonically under the relevant physiological conditions. When considering mix- 
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tures of purely repulsive crowders of different sizes the effects of crowding are additive, and 
the most stabilizing composition will consist completely of the smallest crowder. Therefore 
any maximum in stabilization as a function of the ratio of various components at a fixed 
total packing fraction implies the existence of attractive interactions with at least one of the 



crowders 



32 ] . 
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